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HIGHER  ORDER  THEORY  FOR  LONG-TERM  BEHAVIOR  OF 
EARTH  AND  LUNAR  ORBITERS 


BERNARD  KAUFMAN  and  ROBERT  DASENBROCK 

Operations  Research  Analysts 
Naval  Research  Laboratory ,  Washington ,  D.C. 


Abstract:  A  semianalytical  method  canable  of  analyzing  both  lunar  and 
earth  orbiters  is  presented.  Primary  attention  is  focused  on  predicting  the  evolu¬ 
tion  of  the  orbit  as  affected  by  third-body  perturbations  and  those  of  the  rotat¬ 
ing  primary.  The  singly  averaged  (literal)  equations  of  motion  are  expanded  by 
machine  to  a  high  cider  in  the  parallax  factor  and  the  mean  motion  ratio.  These 
equations  are  numerically  integrated  to  yield  the  orbital  evolution  and  stability 
for  a  wide  range  of  initial  conditions.  In  addition,  a  purely  analytical  method  is 
introduced  to  yield  the  orbital  lifetimes  for  a  special  class  of  orbits. 


INTRODUCTION 

The  design  of  any  mission  to  place  a  satellite  in  orbit  about  a  body  is  a  complex  and 
time  consuming  investigation  of  the  types  of  orbits  that  would  meet  stated  scientific 
objectives.  Orbiters  in  the  earth-moon  system  are  influenced  by  a  very  complex  dynamic 
field.  This  system  is  a  unique  one  because  of  the  relative  masses  of  the  two  bodies.  A 
satellite  in  orbit  about  either  of  these  bodies  undergoes  interacting  perturbation^  that  are 
difficult  to  model.  If  periapsis  is  low  enough,  then  the  asphericity  of  the  central  body 
and  atmospheric  drag  become  important;  if  apoapsis  is  high  enough,  perturbations  caused 
by  the  third  body  will  dominate.  Coupling  between  the  gravity  field  and  the  third-body 
perturbations  can  cause  complicated  behavior  in  certain  characteristics  of  the  orbit.  Further, 
the  effects  of  the  sun  as  a  perturbing  force  cannot  be  neglected. 

Because  of  this  complex  dynamic  field,  the  mission  analysis  phase  of  the  study  must, 
and  will  of  necessity,  include  a  time  history  of  many  orbits  in  order  to  gain  the  maximum 
scientific  data  from  the  final  orbit  chosen.  Such  a  history  should  include  not  only  life¬ 
time  predictions  but  also  a  reasonably  accurate  history  of  all  of  the  orbital  parameters 
for  a  variety  of  initial  conditions.  Because  such  a  variety  of  conditions  will  be  used,  it 
also  becomes  essential  that  the  model  chosen  to  produce  the  time  history  be  not  only 
accurate  but  very  fast.  Any  good  n-body  precision  integration  program  of  the  Cowell 
or  Encke  type  is  capable  of  meeting  the  first  of  these  criteria,  but  certainly  not  the  second. 

To  speed  up  the  computations  involved,  the  standard  approximation  used  has  been 
to  doubly  average  the  disturbing  function  due  to  the  presence  of  a  third  body.  The  dis¬ 
turbing  function  is  first  averaged  over  one  orbit  of  the  satellite,  and  then  over  one  revolu¬ 
tion  of  the  central  body  about  the  disturbing  body.  This  process  eliminates  all  short  and 
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medium-period  terms,  leaving  only  the  long-period  perturbations  for  consideration.  How¬ 
ever,  experience  has  shown  that  this  model  has  limited  use  for  the  earth-moon  system. 

Also,  for  a  planet  such  as  Mars  where  coupling  between  oblateness  and  third  body  can 
be  strong,  the  doubly  averaged  system  sometimes  fails.  A  good  discussion  of  this  model 
can  be  found  in  Ref.  1. 

The  singly  averaged  equations  of  motion  have  been  shown  in  Ref.  2  to  be  highly 
accurate  for  orbiters  of  Mars.  Here  only  the  short-period  terms  are  averaged  out  and  the 
equations  of  motion  retain  the  medium-  and  long-period  terms.  Therefore  the  model  is 
also  valid  for  both  near  and  far  orbiters  of  the  earth  and  moon.  However,  as  shown  in 
that  reference,  the  expansion  of  the  third-body  disturbing  function  is  essentially  in  terms 
of  the  parallax  factor  a/r'.  For  a  100,000-km-high  earth  orbiter,  this  parallax  factor  is 
about  0.25  for  the  moon  as  the  disturbing  body.  The  expansion  in  Ref.  2  was  truncated 
to  retain  only  terms  of  second  order  in  the  parallax  factor,  but  for  high  earth  or  lunar 
orbiters  this  is  not  sufficient.  The  expansion  must  be  to  at  least  fourth  order,  and  for  high 
orbiters  (100,000  km)  should  be  carried  even  further. 

When  the  third-body  terms  are  averaged,  the  assumption  is  sometimes  made  that  the 
disturbing  body  does  not  move  significantly  over  one  orbit  of  the  satellite.  However,  for 
high  orbits  this  assumption  is  clearly  violated  in  the  case  of  the  earth-moon  system.  Thus, 
in  carrying  out  the  expansion,  a  time  rate  of  change  for  all  terms  containing  the  third-body 
position  must  be  included.  This  yields  a  further  expansion  of  the  disturbing  function  in 
terms  of  the  mean  motion  ratio  n'/n  —  the  ratio  of  the  mean  motion  of  the  disturbing 
body  to  that  of  the  satellite. 

To  carry  out  the  expansions  in  terms  of  the  parallax  factor  and  the  mean  motion 
ratio,  and  then  to  average  the  equations  of  motion  over  one  orbit,  requires  an  excessive 
amount  of  algebra  for  the  higher  order  terms.  It  is  almost  impossible  to  proceed  beyond 
fourth  ord  ~  in  parallax  using  hand  techniques.  To  aid  in  the  algebraic  computations,  a 
general  algebraic  manipulation  routine  was  developed  which  is  now  operational  on  a  large- 
scale  computer.  This  program  was  used  to  compute  the  averaged  (literal)  equations  of 
motion  to  eighth  order  in  the  parallax  factor  and  to  second  order  in  the  mean  motion 
ratio  with  the  corresponding  cross  terms  up  to  and  including  fifth  order  in  parallax.  The 
entire  expansion  required  approximately  2  min  of  central  processor  time  on  a  CDC  3800 
computer.  The  output  is  FORTRAN  compatible,  card  punched,  and  directly  insertable 
into  the  variation-of-parameters  program  with  no  human  interaction. 

In  addition  to  the  third-body  effects,  the  gravity  harmonics  of  the  rotating  primary 
must  be  considered.  For  the  moon,  these  equations  may  be  averaged  over  the  orbital 
period  since  the  moon  rotates  slowly.  However,  for  the  more  rapidly  rotating  earth,  this 
analysis  is  invalid  because  the  orbital  mean  motion  may  be  nearly  commensurate  with 
the  rotation  of  the  primary.  To  avoid  this  problem,  Gaussian  quadratures  are  used,  i.e., 
the  equations  of  motion  are  numerically  averaged  from  one-half  orbit  behind  to  one-half 
orbit  ahead  of  the  present  position  of  the  satellite.  The  Greenwich  mean  sidereal  time 
enters  explicitly  in  this  analysis.  These  averaged  rates  are  then  used  in  the  total  variations 
of  the  elements.  At  present,  a  full  7X7  and  4X4  field  is  used  for  the  earth  and  moon, 
respectively.  When  the  tesseral  harmonics  are  not  required,  only  terms  containing  J2, 

</22,  </3,  and  J4  are  used,  and  the  variational  equations  are  calculated  explicitly  without 
using  quadratures. 
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Finally,  a  model  is  also  included  for  atmospheric  drag  and  its  effects  on  the  semi¬ 
major  axis  and  eccentricity.  This  model  assumes  a  nonrotating  atmosphere,  no  lift  forces, 
and  a  drag  force  acting  as  a  negative  tangential  component.  The  density  is  calculated  only 
as  a  function,  of  altitude,  and  the  accelerations  are  averaged  over  one  revolution  of  the 
satellite. 

In  addition  to  this  semianalytic  approach,  a  purely  analytical  model  is  also  developed. 
This  model  concentrates  on  the  long-period  terms  and  applies  to  both  lunar  and  high 
earth  orbiters.  To  derive  the  long-period  equations  of  motion,  the  medium-period  terms 
must  be  removed.  This  can  be  done  provided  that  the  earth  and  moon  move  much  faster 
on  their  respective  paths  than  the  line  of  apsides  of  the  perturbed  orbit  under  question. 
This  condition  is  satisfied  for  all  lunar  orbiters  for  which  the  height  does  not  exceed 
approximately  four  lunar  radii.  For  the  low  lunar  orbiters,  the  fastest  periodic  term  in 
the  equations  of  motion  involves  the  earth  angle,  i.e.,  one  revolution  per  month,  whereas 
the  nodal  angle  and  perilune  angle  rotate  with  periods  of  not  less  than  6  months.  Thus, 
these  medium-period  effects  involving  the  earth  angle  can  be  removed  by  a  von  Zeipel 
transformation. 

For  earth  orbiters,  the  analysis  is  more  complicated.  The  line  of  apsides  rotates 
with  an  angular  velocity  of  about  8°  per  day  for  a  low  orbiter,  whereas  the  sun  and  moon 
move  along  their  respective  paths  with  angular  rates  of  1°  and  13°  per  day,  respectively. 
Thus  these  medium-period  terms  cannot  be  removed  for  the  low  orbiters;  however,  the 
analysis  is  valid  for  the  higher  orbits  and  useful  results  can  be  obtained. 

In  the  absence  of  oblateness,  the  solution  to  the  equations  of  motion  can  be  expressed 
in  closed  form  using  elliptic  integrals.  However,  in  certain  special  cases  involving  initially- 
near-circular  orbits,  the  solution  involves  only  the  elementary  functions.  These  special 
cases  are  extended  to  include  oblateness  by  introducing  a  series  expansion.  The  solutions 
are  accurate  and  yield  results  applicable  to  initially-near-circular  orbits.  These  solutions 
yield  tho  long-period  time  history  of  eccentricity  and  pericenter  position  and  give  the 
orbital  lifetime  for  unstable  orbits.  This  procedure  works  well  for  high  lunar  orbits  and 
can  be  extended  to  high  earth  orbiters. 


MACHINE-AUTOMATED  ALGEBRA 

The  use  of  machine-automated  algebra  in  celestial  mechanics  has  become  increasingly 
popular  in  recent  years.  Due  to  the  high  probability  of  error  which  is  introduced  when 
hand  methods  are  used,  it  was  decided  to  develop  the  equations  of  motion  to  a  high 
order  by  computer.  To  this  end  it  was  decided  to  construct  an  algebraic  manipulation 
program  compatible  with  the  CDC  3800  computer  presently  in  use  at  the  Naval  Research 
Laboratory.  This  program  will  manipulate  the  otherwise  involved  literal  Poisson  series 
occurring  in  classical  perturbation  theory.  The  computerized  operations  include  the 
simplification,  ordering,  negation,  addition,  subtraction,  multiplication,  differentiation, 
and  integration  of  the  trigonometric  series  occurring  in  the  theory.  Other  more  specialized 
routines  include  a  binomial  and  taylor  series  expansion.  The  program  is  written  in  FORTRAN 
and  can  be  used  on  any  machine  possessing  a  FORTRAN  compiler,  with  little  or  no 
modifications. 

The  equations  of  motion  (to  be  described  later)  were  developed  entirely  by  computer. 
These  literal  equations  were  automatically  card  punched  in  FORTRAN-compatible  form 
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and  inserted  directly  into  the  variation  of  parameters  program  with  no  human  interaction. 
Literally  thousands  of  terms  were  involved  in  these  expansions,  with  a  savings  in  time  of 
many  months  and  possibly  years.  The  method  has  the  added  advantage  in  that  trivial 
algebraic  and  keypunching  errors  are  eliminated  as  possible  errors  in  the  analysis.  Remain¬ 
ing  errors  can  be  traced  to  those  of  concept  and  the  programming  of  the  literal  expansions. 
A  detailed  description  of  the  algebraic  manipulation  program  may  be  found  in  Ref.  3. 


VARIATION  OF  PARAMETERS 

The  equations  for  the  variation  of  the  keplerian  elements  are  well  known  and  are 
developed  in  any  good  textbook  on  celestial  mechanics.  Therefore  the  equations  will  be 
listed  here,  without  derivation,  in  the  two  forms  that  were  used  in  the  computer  program. 


Lagrange’s  Planetary  Equations 

da  _  _2_  3F 
dt  m  3Af 


*\ 


m  (1--«2)1/2  (a  -  *2)1/2  M  _  M\ 

dt  ena 2  V  '  3M  9w/ 


dft  _ 


_ 1 _ 3 F 

dt  na2(  1  —  e2)1/2  sin  i  3i 

di__ _ esc  i  {  3 F 

dt  /ja2(l — e2)1/2  \3^ 


cos  ( 


dtp  _  (1  -e2)1/2  3F  _  cos  i  3 F 

dt  ena2  de  na2(  1  —  e2)1^  sin  i  3 i 

dM  m  _  _2  3 F  _  (1  -  e2)  3 F 
dt  n  na  5a  nea 2  5e’ 


(1) 


Gauss’  Form  of  the  Equations 


da  _  2 

dt  ~n(  1  —  e2)i/2 


(r  e  sin  f  +  s) 


de 

It 


_  (1  ~  e2)1/2  sin  f  R  (1  —  e2)1/2  /fl2(l  —  e2)  —  r2\g 
na  ena  2  \  r  / 


dSl  _  r  sin(to  +  f)  w 
(ft  na2(l  —  e2)17S”  sin  i 

dj  ,  +  W 

dt  na2(l  —  e2)1'2 


> 


(2) 


the  disturbing  force  A  is  defined  as 


A  =  RUf  +  SUfl  +  JFUA  (3) 

and  A  is  decomposed  into  components  in  the  radial  (JR),  transverse  (S),  and  orbit  plane 
normal  (W)  directions  where 


UA  =  Ur  X  Ufl. 


THE  THIRD-BODY  DISTURBING  FUNCTION 


The  acceleration  experienced  by  a  satellite  under  the  influence  of  a  point  mass  third 
body  is 


(4) 


where  u'  is  the  gravitational  coefficient  of  the  third  body  and  if,  denotes  the  acceleration 
vector  in  inertial  space.  The  position  vectors  to  the  third  body  and  the  satellite  are  r'  and 
r,  respectively. 


z 


Fig.  1  —  Planet-centered  satellite  geometry 
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Equation  (4)  may  be  expressed  as 


where 


Equation  (5)  can  be  written  in  the  form 

r  =  VF(r.r') 


(5) 


(6) 


where 


and 


o  r  r 

cos  S  =  — r- . 
it 

Now  introduce  the  eccentric  anomaly  E  directly  into  Eq.  (7),  i.e.,  let 

^cos  S  =  8A(cos  E  —  e)  +  SB(1  —  e2)112  sin  E 
and  2 

=  82(1  —  e  cos  E) 2, 

where 


(7) 


8  =  -p  is  the  parallax  factor, 

A  =  P  •  ?', 

B  -  Q  •?',  where  P  and  Q  are  as  shown  in  Fig.  1,  and 
r'  is  the  unit  vector  to  the  third  body. 

The  quantity  F(r',r)  was  expanded  directly  to  the  eighth  order  in  8  by  using  the 
algebraic  manipulation  program  previously  mentioned.  Due  to  the  direct  nature  of  the 
expansion,  the  explicit  expressions  for  the  Legrendre  polynomials  were  not  required. 
The  coefficients  of  each  factor  8n  (where  2  <  n  <  8)  were  automatically  collected,  and 
the  disturbing  function  is  then  obtained  in  the  form 

8 

F(r,r’)=7  y  ZnFn(A,B,e,E) 

n-2 


(8) 
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where  each  Fn(A,B,e,E )  is  of  the  form 

n  n 

Fn(A,B,e,E)  -  Cji,tm  AJ&e*  cos  mE  +  Sjktm  Ate*e8  sin  mE  (9) 
mm  0  m*l 

and  both  CJftim  and  S^Sm  are  obtained  as  rational  integer  coefficients.  The  next  step  is 
to  average  the  disturbing  function  over  one  orbital  period,  i.e., 

Fn(A,B,e)  =  **»<*«  58  F„(l-ecosE)dE.  (10) 

The  equations  of  motion  as  determined  by  the  Lagrangian  planetary  equations  listed  in 
Eq.  (1)  are 


i,  =  ffa).  (11a) 

Substituting  Fn  and  its  derivatives  we  then  have  the  averaged  equations  of  motion  given 
by 


x,  *  fi(A,B,y,e)  (lib) 

where 

xi  is  the  Keplerian  state  vector, 


with 


In  the  first-order  variation  of  parameters  analysis  it  may  be  assumed  that  the  orbital 
elements  a,  e,  i,  SI,  and  w  are  constant  as  the  equations  of  motion  are  averaged.  However, 
the  third-body  position  cannot  be  held  constant. 

Therefore  A,  B,  and  y  are  expanded  about  their  nominal  values,  i.e., 

*  -  ♦#(£)<"  -  ■ "«» 4  H(£)2<"  -  M<>>2 + •  •  • 


fTi  =  R  r' 

bs*0!  ’r 


(12) 
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with  similar  expressions  for  B  and  7.  O'  is  the  third-body  central  angle  and 

M  -I  p 

dd'  ~  n'  P  ‘  dt 


where  dr' /dt  is  given  along  with  r'  by  an  analytical  ephemeris.  The  additional  terms  for 
Gq.  (11a)  representing  the  motion  of  the  disturbing  body  then  become 

5  ‘  bA  90'\  n)(M  dB  b0'(  "  YM  Mq) 


+  lly  +  hi8her  order  terms*  (13) 

Averaging  Eq.  (13)  as  above  and  then  adding  Eq.  (lib)  yields  the  total  averaged  equations 
of  motion. 


Since  the  disturbing  function  is  time  dependent,  the  time  rate  of  change  of  the  semi¬ 
major  axis  has  a  nonzero  average.  Although  no  secular  change  results,  small  but  not  in¬ 
significant  twice-monthly  and  twice-yearly  fluctuations  in  the  semimajor  axis  do  appear. 
The  final  form  of  the  equations  of  motion  is  obtained  in  the  parallax  and  mean  motion 
ratio  as 


*« 


(14) 


where  the  fk 0  are  the  functions  derived  from  the  averaged  disturbing  function  Fn,  and 
hi  is  the  additional  term  due  to  the  motion  of  the  disturbing  body.  The  derivatives  of 
Fn  necessary  for  fk0  are  listed  in  Appendix  A,  and  the  derivatives  required  for  fkl  are 
listed  in  Appendix  B. 


For  a  very  high  earth  or  lunar  orbiter  it  is  not  enough  to  assume  that  the  orbital 
elements  a ,  e,  i,  SI,  and  a>  remain  constant  during  the  averaging  process,  as  was  done  here. 
For  these  high  orbits  it  is  necessary  to  include  the  coupling  between  the  short-period  fluc¬ 
tuations  in  the  elements  with  the  short-periodic  part  of  the  disturbing  function.  This  is 
especially  important  for  h.'gh  lunar  orbiters.  These  additional  expansions  are  being  carried 
out  but  have  not  yet  been  implemented  in  the  variation-of-parameters  program. 


GRAVITATIONAL  FIELD  ANALYSIS 


For  a  satellite  of  the  moon,  the  disturbing  function  could  be  averaged  analytically; 
for  most  other  planets  such  is  not  the  case.  Representing  the  gravitational  field  of  the 
body  by  the  standard  expansion  in  spherical  harmonics,  we  have 


00  n 

^  ^*(sin  <j>)  (C„m  cos  mX  +  Snm  sin  mX) 

n-i  mm 0 


(15) 
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Then  in  the  Lagrangian  equations  we  would  normally  substitute 

F=-f  Fdt 

7  Jo 

for  the  disturbing  function  where,  more  explicitly, 

F  =  ~  |  (1  +  e  cos  f )n_1  P™(  sin  <(>)  e/m(X_0)  <p  df  (16) 

with  6  -  6t  (0  is  the  rotation  rate  of  the  primary). 

The  integral  in  the  above  equation  is  evaluated  with  all  quantities  held  constant, 
except  f.  For  a  satellite  of  the  moon,  the  change  in  6  over  one  orbit  is  negligible  to  the 
first  approximation,  but  for  other  planets,  this  change  can  be  significant.  It  is  for  this 
reason  that  the  disturbing  accelerations  caused  by  a  nonspherical  central  body  must  be 
treated  differently.  Here  the  equations  of  motion  are  averaged  numerically  and  these 
averaged  equations  are  then  integrated  to  obtain  the  variation  in  the  orbital  elements. 

This  treatment  was  suggested  by  Uphoff  in  Ref.  4.  He  defines  the  averaged  Lagrangian 
equations  of  motion  in  the  following  manner:  Let  x,-  be  the  standard  variational  equation 
of  any  element.  The  averaged  Lagrangian  equation  is  then 

i, -y  5  J% (17) 

But 

df  _  (fxP)1/2 . 
dt 

therefore 

-  nP2 

arVJp 

Substituting  the  Gaussian  form  of  the  variational  equations  into  Eq.  (18)  the  following  is 
obtained: 


r2  ’ 


f 


ii(f)df 


(1  +  e  cos  f  )2 


(18) 


-r  o2(l  —  e2)  (f*  e  sin  f  R(f)  +  (1  +  e  cos  f)S(f) 
n  i{ia  J/x  (1  +  e  cos  f )2 


t  _  (f2  /  sin  f  +  2cos  f  +  e(l  +  cos2  f) 

2 irji  \(1  +  e  cos  f )2  (1  +  e  cos  f)3 

nP2  cos(oj  +  f)W(f) 

*  (l  +  e  cos  03 


S(f)]df> 


(19) 
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-r  =  np2  ff2  /(2  +  e  cos  f)  sin  f  _  cos  f 
2rr/i  \  e(l  +  e  cos  f)3  e(l  +  e  cos  f)2 

_  *■.<-->  */);<"'.  gff)W 

(1  +  e  cos  /)3  sin  i  f 


S2  = 


«P2 

2it(x 


— «W(/) 
(1  +  e  cos /)■*  sin  i 


>  (19) 


M  =  n  +  nt>2{1  ~e2)1/2  I*2  fVC0S  £  _  2  \ 

2irju  |_\  e  (1  +  c  cos  f)J 

x  — ML-  _  (LLl£21^iir  S(/)l  dl 

(1  +  c  cos  f)2  e(l  +  e  cos  f)3  J  ^ 

The  evaluation  of  the  definite  integrals  in  these  equations  is  done  using  Gaussian  quadratures 
with  24  points. 

The  computation  of  the  disturbing  accelerations  is  accomplished  using  a  technique 
devised  by  DeWitt  in  Ref.  5.  This  method  was  programmed  to  include  any  order  of  the 
zonals  and  tesserals  desired.  At  present  a  full  7X7  and  4X4  field  is  used  for  the  earth 
and  moon,  respectively.  The  accelerations  are  evaluated  in  cartesian  coordinates  and  trans¬ 
formed  to  components  in  the  radial  (J?),  transverse  (S),  and  orbit  plane  normal  (IV)  direc¬ 
tions  by  means  of  the  rotation  transformations 

tr&w  =  Tz(w  +  f)  r*(0  Tz(&)- 
The  transformed  accelerations  are  then  used  d'rectly  in  Eq.  (19). 

For  central  bodies  other  than  the  earth  or  moon,  only  terms  containing  </2>  <J?2.  J3, 
and  JA  are  used,  and  the  variational  equations  are  calculated  explicitly  without  going 
through  the  above  averaging  and  quadratures.  The  variational  equations  will  simply  be 
listed  here.  The  equations  in  «/2  and  j£  are  derived  in  detail  in  Ref.  6,  with  additional 
terms  from  Ref.  11,  and  those  in  «/3  and  </4  were  taken  from  Ref.  7.  This  option  is  also 
used  for  high-speed  analysis  of  earth  and  lunar  orbiters  where  a  full  gravity  field  is  not 
required. 


V  (20) 

45 nJ%R*  /14  .  „  \  .  0  .  .  „ 

- - —  e(l  —  e*)  ( —  —  sin^  1)  sm^  »  sm  2w 

32 P*  \15  / 


Equations  in  J2,  J2 
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3 nJ2Re  cos  i  9 nJ§R,  cos  i 


2P2 


2  lve 
4i>4 


*  f*3  _  5 
"  2  3 


sin2  i 


A  .  5  .  o  \  .  «2cos  2<o  /7  „  .  0 

e  le  +  24 s,n  ‘)  *— — ir8'”  ') 


+  (1  -c2)1/2  (l  -- 

/di\ 

_  45  <Re  o/14 

\dt)j 

J2 

64  p4  116  \15 

/dw\ 

JnJ2R?  (  6 

\dr  /, 
«•> 

2P2  \  2 

3nJ%R* 
16 P4 


48  —  103  sin2  i  +  -^7^  sin4  i 
4 


+  e2^7  —  ^  sin2  i  —  ~  sin4  ij  —  cos  2w  ^7  —  —  sin2  i^sin2  i 


—  e2(l  — -y  sin2  i  +  sin4 
+  (1  —  e 2)1/2  (24  —  66  sin2  i  +  45  sin4  i) 


(20) 


where 

Re  is  the  equatorial  radius  of  the  planet 

P  =  a(l  —  e2)  is  the  semilatus  rectum  and  the  subscript  J2  means  oblateness  and  J2 
terms  only. 


Equations  in  </3 


3  nR$Ja 

8  p3 


(1  —  e2)  (5  cos2  i  —  1)  cos  w  sin  i 


3  nRfJo 

- — -  e(15  cos2  i  —  11)  sin  w  cot  i 

3  P2 


> 


3nRgJ2 

- - —  e(5  cos2  i  —  1)  cos  w  cos  i 

8 


(21) 
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m  ■ 


_  3 nRp3  (i  +  4^2) 


8P3  « 


(5  cos2  i  —  1)  sin 


.  .  /dS2\ 

Wsm‘ 


cos  t 


(21) 


where 


m  +2ae§) 

/da\  }dth  Wy, 

\dt/r  (1  e2) 


Equations  in  J4 

(f),  '2Pt“‘(f) 


2P  tan  i 
J4  "“'J4 

15  nR*J, 


(-)  « 

\dt  A  ~  32P4 


<?"4 


32P4 
/<K2\  lbnRgJ4 


e(l  —  e2)  (7  cos2  i  —  1)  sin  2u  sin2  i 


cos  i  {2(7  cos2  i  —  3) 


+  e2[7  cos2  *  —  1  +  4  sin2co(7  cos2  i  —  4)]} 
di  \  15  nRgJ4 


/di\ 

UA 


64P4 
15nfl4</4 


16P4 


e2(7  cos2  i  —  1)  sin  2w  sin  2i 

f 

8  —  28  sin2  i  +  21  sin4  i 


—  (7  cos2  i  —  1)  sin2w  sin2i  +  e2 
in2co^6  —  35  sin2  i  +  sin4  ijj  j- 


+  sin 


(22) 


>  (23) 


63 

6  —  14  sin2  i  +  sin4  i 


where  a  definition  identical  to  Eq.  (22)  holds  with  the  proper  change  in  the  subscripts. 
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DRAG 

For  the  variations  due  to  the  presence  of  an  atmosphere,  a  model  has  been  assumed 
in  which  there  are  no  lift  forces  present,  the  drag  force  acts  as  a  negative  tangential  com¬ 
ponent,  and  the  atmosphere  is  nonrotating.  This  results  in  variations  only  in  the  semi¬ 
major  axis  and  the  eccentricity.  The  variational  equations  are  then  averaged  over  a  single 
orbit  using  Gaussian  quadratures.  The  density  is  taken  from  several  models  and  is  calcu¬ 
lated  as  a  function  of  altitude.  These  models  can  be  found  in  Ref.  8  and  9  for  Mars  and 
Venus,  respectively.  For  the  earth,  no  atmosphere  was  used  in  the  present  program, 
although  it  would  be  very  easy  to  incorporate  a  model  similar  to  the  type  used  for  Mars 
and  Venus. 

The  detailed  derivation  may  be  found  in  Ref.  2  with  only  the  results  being  listed  here: 

-  CDAa2(l  —  e2)3/2  rv  Dyz 

d  =  -  — - 1 - - —  - - df  (24) 

2tt  mp  J_„  (i  +  e  Cos  f)2 

PV(e  +  co jf)df  (26) 

2jt  m  J_n  (i  +  e  cos  f)2 

where  CD  is  the  aerodynamic  drag  coefficient,  A  is  the  cross-sectional  area,  and  p  is  the 
density. 

It  is  to  be  stressed  that  only  a  nonrotating  atmosphere  has  been  considered;  otherwise 
di/dt  and  dSl/dt  would  be  nonzero.  The  complexities  of  assuming  an  exponential  density 
profile  have  also  been  bypassed  by  averaging  the  effects  over  one  revolution  of  the  satellite. 


SAMPLE  CASES 

The  method  described  above  has  been  programmed  for  a  CDC  3800  computer  in 
double-precision  mode  under  the  program  name  of  POPLAR  (Planetary  Orbiter  Prediction 
and  Lifetime  Analysis  Routine).  Figure  2  shows  a  comparison  between  POPLAR  and  an 
Encke  n-body  numerical  integration  program  for  a  lunar  orbiter.  Most  of  the  slight  dif¬ 
ferences  noticed  here  are  due  to  the  fact,  that  mean  elements  were  not  used  and  the  con¬ 
stants  were  different  between  the  two  programs.  However,  it  is  important  to  note  that 
these  differences  do  not  appear  to  be  growing  with  time.  The  numerical  integration  pro¬ 
gram  took  4.6  min  of  360/95  time,  which  translates  to  about  115  min  of  3800  time. 
POPLAR  took  3.65  min  for  the  same  case  —  a  factor  of  over  30  to  1. 

Figure  3  is  a  plot  of  eccentricity  versus  argument  of  pericenter  for  a  lunar  orbiter. 
Superimposed  on  this  graph  are  contours  of  constant  lifetime.  The  value  of  such  a  plot 
is  that  it  allows  initial  conditions  to  be  selected  that  will  yield  any  given  lifetime.  For 
lunar  orbiters,  the  inclination  does  not  vary  by  more  than  a  degree  or  two  and,  therefore, 
these  curves  may  be  used  as  a  very  good  first  approximation  to  a  nominal  orbit.  The 
data  shown  in  Fig.  3  is  composed  of  thirteen  different  cases  and  took  a  total  of  about 
15  min  worth  of  3800  CPU  time. 
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LUNAR  0RB1TER 


Fig.  2  —  Comparison  of  POPLAR-generated  and  numerically  integrated 
values  of  eccentricity  vs  time 


LUNAR  ORBITER 
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W (DEG) 

Fig.  3  —  Eccentricity  vs  argument  of  perigee  for  a  lunar  orbiter 

Figure  4  is  a  similar  curve  for  an  earth  orbiter  of  75°  initial  inclination  to  the  ecliptic. 
For  earth  orbiters,  the  inclination  does  not  remain  as  constant  as  for  lunar  orbiters.  How¬ 
ever,  the  change  in  inclination  is  only  about  12°;  therefore  approximate  initial  conditions 
may  still  be  obtained  but  must  be  finally  checked  by  numerical  integration.  Total  CPU 
time  for  Fig.  4  was  about  25  min. 
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W  (DEG) 


Fig,  4  —  Eccentricity  vs  argument  of  perigee  for  an  earth  orbiter 
with  i  =  76° 


CO  (DEG) 


Fig.  5  —  Eccentricity  vs  argument  of  perigee  for  an  earth  orbiter 
with  i  =  45° 


Figure  5  is  again  for  an  earth  orbiter  inclined  at  45°  to  the  ecliptic  and  shows  stable 
orbits  as  far  as  the  lifetime  is  concerned .  The  intersections  of  the  curves  near  oo  =  ±  90° 
is  indicative  of  the  fact  that  inclination  does  not  remain  constant,  varying  as  much  as  15° 
for  some  of  the  curves  in  this  figure.  Again,  however,  first  approximations  to  initial  con¬ 
ditions  may  still  be  obtained.  Approximately  20  min  of  CPU  time  was  required  for  Fig.  5. 


16 


KAUFMAN  AND  DASENBROCK 


ANALYTICAL  RESULTS 

The  results  obtained  thus  far  have  been  arrived  at  by  numerically  integrating  the 
medium-  and  long-period  equations  of  motion.  For  a  lunar  orbiter  (and  some  earth 
orbiters)  the  medium-period  terms  may  be  averaged  out,  leaving  the  long-period  equations 
of  motion.  For  a  low  lunar  orbiter  these  equations  may  be  integrated  analytically  away 
from  resonance  by  the  method  of  successive  approximations,  as  in  Ref.  10.  In  the  absence 
of  oblateness  the  solution  may  be  expressed  in  terms  of  elliptic  integrals  (Ref.  1).  However, 
in  the  case  of  the  initially-near-circular  orbit,  the  solution  to  the  latter  problem  may  be 
expressed  exactly  in  terms  of  the  elementary  functions.  J2  is  considered  by  introducing 
a  power  series  expansion  in  terms  of  the  quantity  1  —  (1  —  e2)1/2,  which  converges 
quickly  for  most  values  of  e. 

The  long-period  Hamiltonian  for  a  moderately  high  lunar  orbiter  can  be  written  as 


F=— 

«  L 

-3»i\ 

nea2  [/ 
+  —  (5 

,  -  3~\ 

fs  S. 

2L2  4 

L3G3  \ 

OV 

16 

t2/ 

\  G2 

(26) 


where  L,G,H,  and  g  are  the  usual  Delaunay  variables.  For  a  Hamiltonian  system 


-  _  0F  .  dF 

G  -  — ;  g  —  ~ 
9 g  9G 


8  \  W\ 


(27) 


(28) 


A  relationship  between  sin  2 g  and  L,G,H  must  now  be  determined.  This  is 
accomplished  by  setting 


F(L,G,H,g) \QmL  =  F(L,G,H,g) 


(29) 


and  solving  for  cos  2g: 


cos  2 g  = 


(-a 

(-5) 


H2 

1  —  3— 

-HU 

-3H\ 

Ko 

X  u 

L2 

G3  \ 

G2/ 

, 

G2 

l1 

1  Lt) 

(30) 
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where 


Kq  is  the  ratio  between  the  lunar  oblateness  effect  and  the  terrestrial  gravity  effect. 
Now  let 


«G>  -  sr 

(31) 

m-± 

(-«) 

(32) 

and  expand  the  expression  for  cos  2g  up  to  the  second  power  of  (1  —  G/L)  =  X: 

cos24Vi-*?>i[1 

5^  (1  +  2X  +  3X2  +  ...)J 

5l1  <32/ 

1  KoL3 

2  7~±m (Ki + K*x + K*X) 

(33) 

where 


and 


„  _  15  H2  1 
Kl  =  ~2  ~~  2 

„  _  105  H2  5 
*2  4  L2  4 


525  H2  55 
'  8  L2  24 


cos  2 g  = 


i40  +  AlX  +  A2X2 


A^l-^\KoL3K> 

U2  1 

A' 


(34) 
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f]2  -i 

*^+5*0*^* 

o “  „!„2  fx  -  gW  *  BtX  *  B2&) 112 .  (35) 

Here 

B0  =  l-p-A§ 

H2  (  H^\ 

B1  =  2A0Al-4jj\l-^ 

ij2  u2  (  f/2\ 

In  terms  of  X  the  equation  for  G  becomes 

X  =  (C0  +  qX  +  C2*2)1/2  (36) 


Equation  (37)  gives  a  time  history  for  the  evolution  of  a  lunar  orbit  on  the  boundary 
contour  separating  the  librating  and  circulating  orbits.  The  solution  is  valid  for  any  value 
for  the  semimajor  axis  as  long  as  the  boundary  situation  occurs.  The  condition  for  the 
existence  of  an  equilibrium  point  (and  boundary  line)  is 
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These  results  can  be  summarized  in  Table  1.  The  lifetimes  for  six  different  lunar 
orbiters  are  computed  using  Eq.  (37).  As  a  check  these  same  cases  are  computed  using  the 
program  POPLAR  which  is  described  in  a  previous  section.  As  can  be  seen  from  Table  1, 
the  agreement  is  excellent  considering  the  simplicity  of  the  analysis. 


Table  1 

Comparison  of  Computed  (Eq.  37)  and  POPLAR-Generated 
Lifetimes  for  Six  Lunar  Orbiters 


Ao 

(km) 

*0 

*0 

(deg) 

w0 

(deg) 

SIq 

(deg) 

Lifetime  (yr) 

Eq. (37) 

POPLAR 

5214 

0.1 

90 

40 

0 

0.92 

0.99 

5214 

0.1 

75 

40 

0 

0.98 

1.05 

5214 

0.2 

75 

40 

0 

0.64 

0.68 

6952 

0.1 

90 

40 

0 

0.65 

0.70 

6952 

0.1 

75 

40 

0 

0.70 

0.79 

6952 

0.2 

75 

40 

0 

0.48 

0.55 

SUMMARY 

The  inclusion  of  medium-period  terms  in  the  equations  of  motion,  while  not  being 
new  in  the  theory,  has  not  until  the  present  been  available  in  a  rapid  program  useful  for 
the  earth  and  moon  because  of  the  excessive  algebra  involved  in  carrying  out  the  expan¬ 
sions  in  the  parallax  factor.  The  development  of  an  algebraic  manipulation  routine  has 
not  only  made  this  expansion  possible,  but  has  allowed  the  further  inclusion  of  the  motion 
of  the  disturbing  body  into  the  Lagrangian  equations.  Usually,  the  disturbing  body  has 
been  held  fixed  during  the  averaging  process.  But  for  high  orbiters  this  assumption  is  no 
longer  valid.  The  algebraic  program  carries  out,  automatically,  the  expansion  of  the  dis¬ 
turbing  function,  differentiation  of  the  expansion,  averaging  of  the  equations  of  motion, 
and  then  punches  out  the  resulting  equations  on  computer  cards  in  a  FORTRAN-compatible 
mode.  These  cards  are  then  inserted  directly  into  the  program. 

Gravitational  harmonics  for  the  earth  and  moon  are  also  included  by  averaging  the 
variations  in  the  elements  by  means  of  Gaussian  quadratures  over  one  orbit  of  the  satellite. 
These  averaged  variations  are  then  included  in  the  total  variation  of  the  elements  which 
are  numerically  integrated  to  yield  a  time  history  of  the  orbit. 

Drag  effects,  while  not  used  in  the  present  examples  for  the  earth,  may  be  included 
in  the  perturbation  model.  At  present,  this  is  limited  to  a  nonrotating  atmosphere  with 
no  lift  present  and  with  the  drag  force  acting  as  a  negative  tangential  component.  The 
atmospheric  density  p  is  calculated  as  a  function  of  altitude. 

The  program  has  proven  to  be  a  very  fast  and  accurate  one.  In  its  high-speed  mode 
of  calculating  third-body  perturbations  and  oblateness  ( J2 ,  </3,  J4)  only,  it  has  reached 
speeds  of  greater  than  500  to  1  over  numerical  integration.  When  the  full  gravitational 
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harmonics  are  included,  speeds  of  25  to  50  to  1  are  still  attainable.  Such  speeds  and 
accuracy  make  the  program  extremely  useful  in  parametric  studies  and  in  gaining  insights 
into  the  behavior  of  planetary  orbiters. 

In  the  theoretical  development,  points  of  singularity  exits  at  i  -  0°  and  180°,  and 
e  =  0.  However,  in  the  practical  application  of  the  theory,  one  may  come  arbitarily  close 
to  these  values  without  encountering  problems. 

In  special  cases  it  has  been  shown  that  the  long-period  equations  of  motions  may  be 
solved  analytically.  The  results  yield  lifetimes  of  unstable  orbits  accurate  to  about  10%. 
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Appendix  A 


THE  AVERAGED  THIRD-BODY  DISTURBING  FUNCTION 
AND  ITS  DERIVATIVES 


The  derivatives  listed  below  are  those  required  for  *n  Eq.  (14). 

3F  =  3F  3A  +  9F  9B 
9co  dA  9 to  9B  9co 
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9F 

9  i  dA  9  i  +  9 B  9  i 
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9F  order  •=• 
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disturbing  function  for  ORDER 

F  -  3#£**2*A**2 

-3#E#*2*B**2 
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3*A**2 
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2 

-3#E**i 

2 


BF/BA  FOR  ORDER  2 
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DISTURBING  FUNCTION  FOR  ORDER  3 


F  =  -25*E**3*A**3 
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&F/3B  FOR  ORDER  3 
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DISTURBING  FUNCTION  FOR  ORDER  4 
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8 
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DISTURBING  FUNCTION  FOR  ORDER  5 

F  =  -441*E**5*A#*5 
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DISTURBING  FUNCTION  FOR  ORDER  6 
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5F/BE  FOR  ORDER  6 
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disturbing  function  for  order  7 

F  =  -3861*E**7*A**7 

32 
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-135135*E**7*A**3*B**4 
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-31185*E**7*A*#3*B**2 
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64 
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2048 
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184275#E**5*A**1*B**2 

2048 

11025*E**5*A**1 

256 

-135135*E**3*A**7 

256 

-1756755*E#*3*A**5*B**2 

2048 

18814V5*E**3*A**5 
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-135135*E**3*A**3*B*#4 
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738045*E**3*A**3*B**2 
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64 
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31185*E#*1*A**1*B**4 
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12b 

-99225*E**6*A**3 

256 

945945*E**6*A**1*B**6 

2048 
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~3690225*E**4*A**3*B**2 
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31185*A**5 
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~405405*A**3*6**4 
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256 

-1216215*E**1*A**2*B**4 

2048 

93555*E**1*A*#2*B**2 

128 

-25515*E**1*A*#2 

128 

-135135#E**1*B**6 

2048 

31185#E**1*B**4 

256 

-8505*E**1*B**2 

128 

315*E**1 

32 

8F/5B  FOR  ORDtK  7 

81081*E**7*A**5*B**1 

64 

-135135*E**7*A*#3*B**3 

64 

-31185#E**7#A**3*B**1 

64 

405405*E**7*A*#1*B##5 

1024 

93555#E**7*A**1*B##3 

512 
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42525*E**7*A**1*B**1 

1024 

27027#E**5*A**5*B**1 

32 

1756755*E**5*A*#3*B**3 

512 

-738045*E**5*A**3*B**1 

512 

*■12162 15*E**5*A**1*B**5 
1024 

31185*E**5*A*#1*B**3 

256 

184275*E**5*A**1*B«*1 

1024 

-1756755*E**3*A**5*B**1 

1024 

-*135l35*Ei  *3*A**3*B**3 
256 

738045*E*  *3#A**3*B**1 
512 

1216215*E**3*A**1*B**5 

1024 

-405405*E**3*A**1*B**3 

512 

-2835*E**3*A##1*B*#1 

32 

'•405405*E**1*A##5*B**1 

1024 

-405405*E**1*A**3*B**3 

512 

31iy5#E**l#A**3*B**l 

64 

-405405*E**1*A**1*B**5 

1024 

31185*E»*1*A##1*B**3 

64 

-8505*E**1*A##1*B**1 

64 
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DISTURBING  FUNCTION  FOR  ORDER  8 


F  =  32175*E**8*A*#8 

128 

-225225*E**8*A**6*8**2 

128 

-45045*E**8*A#*6 

128 

1126125*E**8*A*#4*B**4 

512 

225225*E**8*A*#4*B**2 

256 

86625*E**8#A*#4 

512 

— 11261 25*E**8*A**2#B#*6 
2048 

-  675675*E**8*A**2*B**4 

2048 

—  259875*E**8*A**2*B**2 

2048 

-55125*E**8*A**2 

2048 

225225,»E*'*8*B*,*8 

16384 

45Q45*E#*8#B**6 

4096 

51975*E**8#B##4 

8192 

r  11025*E**8*B**2 

4096 

11025*E*#8 

16384 

225225*E**6*A**8 

128 
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-  675675*E**6*A**6*B**2 
256 

-795795*E**6*A**6 

256 

-5630625*E**6*A**4*B**4 

2048 

3828825*E**6*A**4*B**2 

1024 

3378375*E**6*A**4 

2048 

6531525*E#*6*A**2*B**6 

4096 

-2027025*E**6*A**2*B**4 

4096 

-3066525*E**6*A**2*B**2 

4096 

-1135575*E**6*A*#2 

4096 

-  225225*E**6*B**8 
4096 

-15Q15*E**6*B**6 

4096 

86625*E**6*B**4 

4096 

77175*E**6*B*#2 

4096 

3675*E**6 

512 

1126125*E**4*A**8 

512 

5630625*t**4#A**6*B**2 

2048 

-8783775*E**4*A**6 

2048 

-8333325*E**4*A**4*B**4 

8192 
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-9594585*E**4*A**4*B**2 

4096 

21008295*£**4*A*#4 

8192 

-6081Q75*E**4*A**2*B**6 

4096 

1936935*E**4*A**2*B**4 

1024 

426195*E**4*A**2*B**2 

4096 

-240345*£**4*A*#2 

512 

675675*E**4*B**8 

8192 

-  225225*E**4*B**6 
4096 

-336105*E**4*B*#4 

8192 

2205#E**4*B**2 

512 

6615*E**4 

512 

1126125*E**2*A**8 

2048 

6531525*E#*2#A**6*B**2 

4096 

-4639635*E**2*A**6 

4096 

6081075#E**2*A##4*B**4 

4096 

-8963955*E**2*A**4*B**2 

4096 

377685*E#*2*A**4 

512 
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1576575*E**2*A**2*B**6 

4096 

~4009005*E**2*A*#2*B**4 

4096 

93555*E**2*A**2*B**2 

128 

“40005*E**2*A**2 

256 

—  225225*E**2*B**8 
4096 

3153l5*E**2*B**6 

4096 

-3465*E**2*B#*4 

512 


-5355*E**2*B**2 

256 

315*E**2 

64 


225225*A**8 

16384 

225225*A*#6*B*#2 

4096 

-15015*A**6 

512 

675675*A**4*B*#4 

8192 

-45045*A**4*B**2 

512 

10395*A**4 

512 

225225*A**2*B*#6 

4096 

-45045*A**2*B*#4 

512 

10395*A#*2*B**2 

256 
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-315*A**2 

64 

225225*B**8 

16384 

-15015*8**6 

512 

10393*B**4 

512 


315*B**2 

64 

35 

128 


BF/3E  FOR  ORDER  8 

32l75*E**7*A*#8 

16 

-225225*E**7*A**6*B**2 

16 

-45045*E**7*A**6 

16 

1126125*E**7*A**4*B**4 

64 

225225*E**7*A**4*B**2 

32 

86625*E**7*A**4 

64 

-1126123*E**7*A**2*B**6 

256 

S’  — 675675*E**7*A**2*B**4 

256 

-  259875*E**7*A**2*B**2 
256 

-55125*E**7*A**2 

256 


NRL  REPORT  7527 


225225*E**7*B**8 

2048 

45045*E**7*B*#6 

512 

51975*E**7*B**4 

1024 

11Q25*E**7*B**2 

512 

1 1025#E**7 
2048 

675675*E**5*A*#8 

64 

-2027025*E**5*A*#6*B**2 

128 

-2387385*E**5*A**6 

128 

-16891876*E**5*A**4*B**4 

1024 

11486475*E**5*A**4*B#*2 

512 

10135125*E**5*A**4 

1024 

19594575*E**5*A#*2*B#*6 

2048 

-6081075*E**5*A**2*B**4 

2048 

-9199575*E**5*A*#2*B**2 

2048 

-3406725*E**5#A*#2 

2048 

-675675*E**5*B**8 

2048 

-45046*E**5*B**6 

2048 
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259875#E**5*B**4 

2048 

231525*E**5*B*#2 

2048 

11025*E**5 

256 

1126125*E**3*A**8 

128 

5630625*E**3*A*#6*B**2 

512 

-8783775*E**3*A**6 

512 

-8333325*E**3*A**4*B**4 

2048 

-9594585*E**3*A**4*B**2 

1024 

21008295*E**3*A**4 

2048 

-6081075*E**3*A**2*B**6 

1024 

1936935*E**3*A**2*B**4 

256 

426195*E#*3*A**2*B**2 

1024 

-240345#E**3*A*#2 

128 

675675*E**3*B**8 

2048 

-225225*E**3*B**6 

1024 

-336105*E**3*B**4 

2048 

2205*E**3*B**2 

128 


NRL  REPORT  7527 


6615*E**3 

128 

1126125*E**1*A**8 

1024 

6531525*E#*1*A**6*B#*2 

2048 

-4639635*E**1*A*#6 

2048 

6081Q75*E**1*A**4*B**4 

2048 

-8963955*E**1*A**4*B**2 

2048 

377685*E**1*A#*4 

256 

157657:>*E*»1*A*#2*B**6 

2048 

-4009005*E**1*A**2*B**4 

2048 

93555*E**1*A**2*B**2 

64 

-40005*E**1*A*#2 

128 

-225225*E**1*B**8 

2048 

3153l5*E**l*B**6 

2048 

-3465*E**1*B*#4 

256 

-5355*E**1*B**2 

128 

315*E#*1 

32 
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SF/3A  for  ORDtR  8 


32175*E**8*A**7 

16 

-675675*E**8*A**5*B**2 

64 

-135135*E**8*A**5 

64 

1126125*E**8*A**3*B**4 

128 

225225*E**8*A**3*B**2 

64 

86625*E**8*A#*3 

128 

-1126123*E**8*A**1*B**6 

1024 

-  675675*E**8*A**1*B**4 

1024 

-  259875*E**8*A**1*B**2 

1024 

-55125*E**8*A**1 

1024 

225225#E**6*A**7 

16 

-2027025*E**6*A**5*B#*2 

128 

-2387385*E**6*A#*5 

128 

-5630625*E**6*A**3*B**4 


3828825*E**6*A**3*B**2 

256 

3378375*E**6*A**3 

512 
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6531525*E**6*A**1*B**6 

2048 

-2027025*E**6*A**1*B**4 

2048 

-3066525*E**6*A**1*B#*2 

2048 


-1135575*E#*6*A*#J 

2048 

1126125*E**4*A*#7 

64 

16891875*E**4*A**5*B*#2 

1024 

“■2635l32i>#E**4*A#*5 

1024 

-8333325*E**4*A**3*B**4 

2048 

-9594585*E**4*A**3*B**2 

1024 

21008295*E**4*A**3 

2048 

-6081075*E**4*A**1*B**6 

2048 

1936935#E**4*A**1*B**4 

512 

426195*E**4*A**1*B**2 

2048 


-240345*E**4*A**1 

256 

1126125*E**2*A*#7 

^  256 

19594575#E**2*A**5#B**2 

2048 

-13918905*E#*2*A*#5 

2048 
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6081075*E**2*A**3*B**4 

1024 

-8963955*E**2*A**3*B**2 

1024 

377685*E**2*A**3 

128 

1576575*E**2*A**1*B**6 

2048 

-4009005*E**2*A**1*B**4 

2048 

93555*E**2*A**1*B**2 

64 

-40005*E**2*A**1 

128 

225225*A**7 

2048 

675675*4**5*8**2 

2048 

-45045*A**5 

256 

675675*A**3*B**4 

2048 

-45045*A**3*B**2 

128 

10395*A**3 

128 

225225*A**1*B**6 

2048 

-45045*A**1*B**4 

r*  256 

10395*A**1*B**2 

128 

-3l5*A**l 

32 
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9F/9B  FOR  ORDtR  8 


-225225*E**8*A**6*B**1 

64 

1126125*E**8*A**4*B**3 

128 

225225*E**8*A**4*B**1 

128 

“3378375*E**8*A*#2*B**5 

1024 

-  675675*E**8*A*#2*B#*3 
512 

-259875*E**8*A**2*B**1 

1024 


225225#E**8#B*#7 

2048 

135135*E**8*B*#5 

2048 

51975#E**8*B**3 

2048 

11025*E**8*B*#1 

2048 

-  675675*E**6*A*#6*B**1 
128 

“5630625*E**6*A*#4*B**3 

512 

3828825*E**6*A*#4*B**1 

512 

19594575*E**6*A**2*B**5 

2048 

-2027025*E**6*A*#2*B**3 

1024 

-3066525*E**6*A*#2*B**1 

2048 
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-225225»E**6*B**7 

512 

-45045*E**6*B**5 

2048 

86625*E**6*B**3 

1024 

77175*E**6*B**1 

2048 

5630625*E**4*A*#6*B**1 

1024 

-8333325*E**4*A*#4*B**3 

2048 

-9594585*E**4*A*#4*B**1 

2048 

“182432 25*E**4#A*#2*B**5 
2048 

1936935*E*#4*A**2*B**3 

256 

426195*E**4*A**2*B**1 

2048 

675675*E**4*B**7 

1024 

-67567S*E**4*B**5 

2048 

-336l05*E**4*B**3 

2048 

2205*E**4*B**1 

256 

6531525*E**2*A**6*B**1 

2048 

6081Q75*E**2*A**4*B#*3 

1024 

-8963955*E**2*A*#4*B**1 

2048 
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4729725*E**2*A**2*B**5 

2048 

-4009005*E**2*A**2*B**3 

1024 

93555*E**2*A**2*B**1 

64 

-225225*E**2*B**7 

512 

945945*E**2*B**5 

2048 

-3465*E**2*B**3 

128 

-5355*E**2*B**1 

128 

225225*A**6*B**1 

2048 

675675*A**4*B*#3 

2048 

“45045*A**4*B**1 

256 

675675*A**2*B**5 

2048 

“45045*A**2*B**3 

128 

10395*A**2*B**1 

128 

225225*8**7 

2048 

“45045*8**5 

256 

10395*8**3 

128 

“315*8**1 

32 
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THE  FOLLOWING  IS  FOK  OKDEK  2 

D1  =  -9*t**2*A**l 

4 

tt*t**l*A**l 

3*A**1 

2 


D2  =  y#E**2*b**l 

4 

-8*E**l*b**l 

-3*13**  1 
2 


D3  =  -3*t**2*A**l 

S*E**1*A**1 

3*A**1 

2 


r 


D4  =  3*E**2*B**1 

2 


3*tS**l 

2 
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D5  =  3#E**l*b*#l 

4 

-3#B**1 


D6  =  -9*E#*3*B**1 

0 

4*E*#2#U**1 

3*t*#l*B**l 

4 


D7  =  3*E**1*A#*1 

4 

-3*A**1 


D8  =  -V*E**3*A*#1 

8 

4#E#*2*A*#1 

3#fc**l*A**l 
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D9  =  y*c:**2*tS**l 

a 

4 


DIO  =  0 


Dll  = 

8 

-4*E**1 

-a 

4 

D12  =  V*t#*2*A**i 

a 

“4*t**l*A*#l 

-3*A**1 

4 
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THt ‘ HOLLOWING  IS  f-OK  OkDtK  3 


D1  =  4i>*E**3*A**2 

4 

-61i>*E**3#b##2 

64 

-21*E**3 

64 

-76*E**2*A**2 

2 

30*E#*2*d**2 

3*E**2 

2 

-7o6*t**l*A**2 

64 

64S*t*#l*ti**2 

64 

3*t**l 

a 

-16*A**2 

2 

3 

2 


D2  =  -6lt>*E**3*A**l*B**l 

32 

60*t**<:*A**l*b**l 


645*E**1#A**1*B**1 

32 
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r 


D3  =  4b*E**3*A**2 

4 

-76*£*#3*b**2 

16 

-21*E**3 

16 

-43*E##2*A*#2 

2 

S>*E**2 

2 

-13:>*t**l*A**2 

16 

7b*E**l*b**2 

16 

3#£**1 

4 

-16*A**2 

2 

3 

2 

D4  =  A##l*d*#l 

« 

7b*b#*i’*A*#l#b**l 

6 


D5  =  -4a*t#*2*A**l*b*#l 

4 

40^t**l'K'A**l*b**l 

ip*A*#l,'b*^l 

2 
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D6  =  16#E*#4*A**1*B**1 

2 

-2i>*t**3*A**l*B**l 

-23l>*E**2*A**i*B**i 

32 

-3#k#*l*A*#l  •*13**1 


D7  =  -4!>*E#*2*A*#2 

8 

133*£#*2#B**2 

32 

9*E#*2 

32 

20*t*#l*A*#2 

-10*£**l*b#*2 

-2*E**1 

15#A**2 

4 

-3 

4 


D9  =  13*E*#4*A**2 

4 

-133#t:**4*B**2 

64 

“21*E##4 

64 
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-2:>*E**3*A*#2 

2 

5#E*#3 *13**2 

3*E**3 

2 

-23i>*E**2*  A**2 
64 

135#E**2*B**2 

64 

3*E**2 

8 

-3*E**1*A**2 

2 

-;>*E**1*B**2 

3*t**l 

2 


D9  =  -13*t**3*A**:L*U*U 

2 

25*E**2#A**1*B**1 

235*E**1*A**1*B**1 

32 

6*A**1*B**1 


DIO  =  -13*L**3*b**l 

2 

25*E**2*B**1 

23Jj*E**l*B**l 

32 


3*b**l 
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Dll  =  -1:>*E**3*A**1 

2 

2i>*E**2*A##l 

2p3*E**1>!,A##1 

32 

6*A**1 

D12  -lb#e.**3*A#*2 

4 

13iJ*E**3*B**2 

64 


21*E##3 

64 

2:>*E*#2*A#*2 

2 

-6#ii**2*ti*#2 

“3*t**2 

2 

26‘>*E**1*A**2 

64 

-136*E**1*B**2 

64 

-3*t**l 

b 

b*A**2 

2 

B#B#*2 

-3 

2 

D13  =  13i>*E**3#B**l 

32 

-10*E**2*ti**l 

-133*E**1*B**1 

32 


10#b**l 
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THE  FOLLOWING  IS  FOR  ORDtK  4 


D1  =  -175*E**4*A**3 

4 

32i>3*E**4*A**l#B*#2 

32 

13;,#E**4*A#*1 

32 

140*E**3*A**3 

-2y4*E^*3*A*#i*B **2 

-18*E*#3*A**1 

b7i>#E**2*A**3 

16 

-4Qyi>*E**2*A**l*B**2 

32 

-16S*E#*2*A*n 

32 

b4#E**l#A*#3 

-42*E**1*A**1*B**2 

-3Q#E**1*A##1 

1 75#A**3 
16 

-10t)*A**l*B**2 

16 

-15*A**1 

4 


D2  =  3255*E**4*A**2*B**1 

32 

-385*E**4*B*#3 

16 


NRL  REPORT  7627 


-135*E*#4*B**1 

32 

-294*E**3*A**2*S**1 

56*E**3*B**3 

18*E**3#B**1 

-409fj*E#*2*A**2*B**l 

32 

245#E#*2*B#*3 

8 

165#E**2#B**1 

32 

“42*E##1#A*#2*B#*1 

-i>6*E**l*B**3 

30*E**1#B#*1 

-10i>*A**2*B**l 

16 

-1G5*B**3 

16 

li>#B#*l 

4 


-3i>*E**4#A**3 

315#E#*4*A**1#B*#2 

8 

75#E#*4#A##1 

8 


7Q*t**3*A**3 
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-30*E**3*A**1 

1Q5*E**2*A*#3 

4 

-525*£**2*A**1*B**2 

16 

-105*E**;i"A**l 

16 

7q*E**1*A**3 

-30#E#*1*A**1 

173*A**3 

16 

-105*A**1*B**2 

16 

-15*A**1 

4 


D4  =  315*£**4*A**2*B«*1 

8 

-105*E**4*B*#3 

16 

-45*E**4*B**1 

16 

-525*E**2*A**2*B**1 

16 

105*E**2*B**3 

8 

-15*E**2*B**1 

16 

-105*A**2*B**1 

16 

-105*B**3 

16 

1&*B**1 

4 
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315*E**3*A**2*B«*1 

4 

“1085#E#*3*B**3 

64 

-255*E**3*B*#1 

64 

-525*E**2*A**2*B**i 

2 

119#E#*2*B**3 

3 

41*E**2*B##1 

2 

~53b5#E**l*A#*2*B**l 

64 

525*E**1*B**3 

64 

135*E**1*B*#1 

16 
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